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Abstract 

All available data on the total cross sections and i?-ratio of e~^e~ — > hadrons 
are compiled from the PPDS (DataGuide), PPDS (ReacData) and HEPDATA 
(Reaction) databases and transformed to a compilation of data on the e~^e~ 
quark quark — > hadrons in the continuum region which can be used in tests of the 
parton model and pQCD calculations. This evaluated data compilation is made 
available in PPDS system and is accessible through the Web. It is shown that cur- 
rent predictions from the parton model and pQCD are well supported by this world 
"continuum" data compilation which can then be used in future refinements of the 
CKs(Q^) as well as aqEDiQ"^) evolution forms. 



Motivation 

Despite the tremendous efforts of experimentalists and phenomenologists devoted to 
the investigation of hadron production in e~^e~ collisions there is no real integrated 
view on the experimental situation even for one of the main observables, the total 
cross section for the reaction e~^e~ hadrons and the ratio 

^ (7(e+e^ hadrons) 
cr(e+e~ fJ^+H')QED 
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Some data are published only as cross sections, a, and others only as R. The situ- 
ation is more complicated as QED radiative effects to the data have been corrected 
for in different ways, and to different degrees, in the various data sets. For example, 
in constructing R experimentalists have used two forms of theoretical calculations, 
one with agEniO), or another taking into account the then currently accepted form 
of evolution of ctQEDiQ"^)- Furthermore, some data were never published in numer- 
ical form and have been read by compilers from graphs. The current experimental 
situation (March 2001) is summarized in Figures 3 and 4. 

This aim of this work is firstly to prepare a comprehenesive compilation of a 
and R transformed, wherever possible, to a unique, and meaningful, style of im- 
plementing the pure QED radiative corrections. The second aim is to construct 
a procedure to extract "continuum data sets" where the direct comparison of the 
parton model with different variants of the pQCD corrections will be conclusive. In 
selecting data for the compilation we use expert assessments from recent reviews 
and works dealing with precision estimation of running aqEDiQ^) and as^Q"^) (see 

i, 0, i, i, 0, 0, m, m, id- 



1 Data selection and normalizations 

1.1 General overview 



The data, which have been published in [^]-|112], were extracted in numerical form 
from the PPDS RD database and from the Reaction database of HEPDATA 
system [114]. The following criteria were then applied to define suitable data sets: 



• All preliminary data are excluded; 

• Data obtained under incomplete kinematical conditions, and not extrapolated 
to the complete kinematic region, are excluded; 

• Data in the form of dense energy scans are omitted if the authors published 
in addition the data averaged over wider energy bins. 

The data sets remaining after the above exclusion criteria have been applied, 
are subdivided into the following six cathegories: 

d 



(Refs. [pO| - |67]). Cross section data corrected by their authors for the con- 



tributions of two-photon exchange diagrams, the initial state bremsstrahlung and 
for the initial state vertex loops, but not corrected for the vacuum polarization con- 
tribution to the running oqed- These "semi-dressed" cross sections correspond 
to the contribution of the diagram shown on Fig. 1. The full j/Z propagator, taking 
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into account all vacuum polarization effects, is denoted on the figure by the bold 
waved line. 



(J2 (Refs. - Cross section data radiatively corrected by their authors 

according to the procedure of Bonneau and Martin |||]. This procedure included 
radiative corrections for the initial state radiation, electronic vertex correction, and 
the correction for the electronic loop in the photon propagator. Thus it partially 
took into account the vacuum polarization effects also. To obtain from these cross 
sections the "semi-dressed" ones (corresponding to the diagram on Fig. 1) we 
rescale cr^^ by a factor 1/(1 — AaQ£;^(s))^. We suggest that all the cross sections 
published before 1978 were radiatively corrected according to this procedure. 




hadrons hadrons 
Fig. 1 Fig. 2 



R^^^ (Ref. [^). Data on R obtained from measurements of the "semi-dressed" 
{(T^) cross-section divided, by their authors, by the point-like muonic cross section 
with fixed oqed = o:qed{0) We rescale these data by the factor {aQEoiO) / ctQEois))'^ / (1— 
AaQ^£,(s))^, to obtain the "bare" i?-parameter described by the diagram shown 
on the Fig. 2. The tree-level j/Z propagator is denoted here by the waved line. 

i?^"'''' (Refs. m - 11). Data on R obtained fr om measurements of the "semi- 
dressed" (cr^*^) cross section divided, by their authors, by the point-like muonic 
cross section, but with the running aQED{s) thus obtaining the "bare" i?-ratio, 
which corresponds to the Fig. 2 diagram again. 

(Jg^ (Refs. ^l| - [^]). LEP I cross section data at the Z peak not corrected 
for initial state radiation and electronic QED vertex loops, and LEP II - III cross 
section data at ^/s > 130 GeV with a cut s' / s = 1 — {y/s/2){l — 0.7225). One can 
interprete ^fs' as an effective mass of the propagator after the initial state radia- 
tion has reduced the e^e~ pair center-of-mass energy. Indeed, the definition of s' 
depends on the assumptions about the initial and final state radiation interference. 
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These data were rescaled to the "semi-dressed" cross section, cr^*^ , defined by th£ 
contribution of the Fig. 1 diagram as follows. First, the theoretical cross sectionsEI 
cr*Jjj and cr^Born ^^^^ calculated using the ZFITTER 6.30 packagaHl ||2^ assuming the 
recommended by PDG ||2| values of the standard model parameters. cr*(jj means here 
the cross section measured with cuts applied in the particular experiment and calcu- 
lated by ZFITTER. cr^Bom denotes theoretical cross section, calculated by ZFITTER 
in the Improved Born Approximation (IB A), that corresponds to the Fig 1 diagram. 
Then the experimental cross section was multiplied by the factor cr^Boml ^^cut- The 
results of this rescaling procedure are enlisted in the Tables 3a, 3b. 

(Jq^. Low energy data (2m^ < Ecm < 2 GeV). The treatment of these data is 
discussed in the following section. 

1.2 Low energy data treatment (2m7r < Ecm < 2 GeV) 

Our treatment of the cross section, cTq'^ , data in low energy range 2771-,^ < y/s < 2 
GeV requires special consideration. 

All the cross sections a{e~^e~ — > hadrons) in this range are obtained via exclusive 
channel summation and therefore give only a lower estimate of the total hadronic 
cross section. Below 1 GeV we summed up the 2ti and 3-7r channels from the ref- 



erences [q^] - [95 1 using a linear interpolation of the individual data sets within 
specific energy regions and combined them to give the total hadron cross sections 
in these regions. The errors were calculated according to these interpolation and 
summation procedures. 

Such an approach works well if all the data are evenly distributed over and 
have comparable errors. Otherwise if in the given ^/s interval there are few points 
of the leading channel with large errors, and this interval is filled more densely by 
the points of less contributing channels, resonance-like false structures in the total 
cross section may arise. 

The possibility of such false structures is likely for the exclusive channel data 
in the range 1.4 < ^/s < 2 GeV, where we summed the contributions from the 
channels > Shadrons, vr+vr^, K^K~, and KsK^. To avoid these false structures 
we have determined the sum of exclusive channels not for all points, where at least 
one channel is measured, but in a fewer number of points, in which the channels 
yielding the major contributions are measured. In other aspects the summation 



Values for cr*„j were calculated using ZFITTER subroutine ZUTHSM with argument settings accord- 
ing to the cuts applied in the LEP I-II-III measurements. Values for cr^orn were calculated by the same 
subroutine but with the special flag switching it to calculate the IBA cross sections 



2 ZFITTER 6.30 Fortran code is available at |http://www.ifh.dc/~riemann/Zfitter/zfitr6_30Tuu 
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procedure was the same as the one used for the data below 1 GeV. The exclusive 
channel data for 1.4 < y/s < 2 GeV are |9C] - [112|- The exclusive sum data in the 
0.81 < y/s < 1.4 GeV region are taken from |112| ]. The data used for the exclusive 
channel summation around (j) resonance (0.997 < \/s < 1.028 GeV) were corrected 
by their authors for the initial state radiation, electronic vertex loops and leptonic 
1-loop insertions insertions into the 7 propagator. We have properly rescaled these 
data points to the Fig 1 diagram. In the remaining 0.81-1.4 GeV data their au- 
thors applied QED-corrections using the Bonneau and Martin |l| prescription and 
all these data points were rescaled as in the cr^'^-case. 



2 Data compilations on a^^ and R^^^^ 

After the data selection and rescaling, described in Section 1, we assemble the com- 
plete data set of the total hadronic cross sections a^'^, normalized to the contribution 
of the Fig. 1 diagram in accordance with the symbolic relation 



—sd _sd I I _sd I I _sd , , _sd , , 

cr = (Tj^ U <T^ U <Tg U CTg U 



f rfbare , , r}bare\ 



3 4 / QED pole, running oiqed 

and the complete data set for the ii-ratios 

ybare _ 0-{e+e' — ^ Hadvons) pig.^L diagram 



R" 



cl^e^e >■ /J'^ l-l )qED pole, running oiqed 



As decribed above, in creating compilations which are uniform in the sense of a 
standardized implementation of pure QED radiative corrections to the "raw" pub- 
lished experimental data, the values of the running aqEois) and AaQ^^(s) were 
used. The procedures for obtaining their numerical values and estimates for their 
uncertainties are described in the Appendix. Both compilations are stored and 
maintained in the PPDS CS database. Each record in the compilation corresponds 
to the record of the original data but contains the rescaled data points. Brief de- 
scriptions of the applied conversion procedures are stored in the special comment in 
each record. The cr and R^"-^'^ compilations obtained by this procedure are shown 
in the Figs. 5 and 6, respectively. 
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3 The continuum regions 



To be able to pick out the region where the parton model with QCD corrections 
can be tested we need to consider the relevant theoretical formulae. 



3.1 Theoretical relations. 

In the parton model, before QCD corrections are applied, the i?-ratio is given by 



(1) 



where 



with 



Kg 
Kg 



elcg + 2eeeqVeVg Rex + {vl + al)vl |xP 

/_2 , _2\— 2 I 1 2 

{v^ + a^)a \x\ 



^ " 16s2c2 {s- M^ + iMzTz) 
V = ^{T^-2Q-f) 

3\/2Gf 2 , 



(2) 



(3) 



Here , 6^ are an effective sin^ cos^ dw defined through renormalized couplings 

vr • aQED{Ml) 



at s = M|: 



The dominant correction term in p originates from t-quark loops in W and Z prop- 
agators which result in an SU{2) violation due to the large mass splitting between 
the h and t quarks fl^ . 

Including QCD loops the expression for R is now (see, e.g. |21], p2| 



R 



3E 



^ \ vr V vr 



(1.9857 - 0.1153iV/) 



a., 



vr 



(6.6368 + 1.2001iV/ + 0.0052A^? + 1.2395(1; Qq 



(5) 
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The coefficients /i and /2 depend only on the quark velocity /3 = (1 — AM"^ / sY^"^ ^ 
where Mq is the flavour production threshold mass of the quark q. Schwinger [0] 
calculated /i in the QED context and parameterised it in a form which is accurate 
enough for our purposes: 



1 + 1(1 



(6) 



4 '2 

Here A'^^ is a number of active quark flavours. The coefficient /2 (which have no 



counterpart in QED) were calculated by Jersak et al |16]. As there is no compact 
analytical expression for /2, we parameterised it as 



/2 = 04/3^ + «3/?a + (1 - 04 - 03)/?, 



(7) 



with 04 = —16, a?j = 17. ^From equations (|T|), (|2|), (||) one can easily estimate 
that even assuming a 100% error for /2, the relative error of R is less than 10"'^ at 

= 10 GeV and less than 0.5% at the Z pole. As we did not used Z pole data in 
our fits, it could not be a source of large theoretical errors. 

In the massless quark limit, /i = /2 = 1- However, at ^/s = 35 GeV for h quarks 
one has (3 = 0.963 whence /i ~ 1.3 and /2 — 1.7. Indeed, such a parameterisation of 
mass effects is valid in QCD only for 0{as) order. Now the correct parameterisation 
is known upto 0{al) order (|21],|2^), but it was not implemented yet into our 
programme. It is likely that it results in large enough discrepancies just in the 
region where ag/vr becomes large. 

The following three-loop parameterisation of as was chosen here (see, e.g., Ill 



47r 



1/?1 , 



1 

A 



-r ^ tan 



\(3oas 
PI 



PI. 
An 



(8) 



with 



Po 
Pi 
P2 

A 



11 - 2iV//3 , 
2(51 - 19iV//3) , 

(2857 - 50337V/ /9 + 325iV| /27) /2 , 
^AP2/Pl-P!/Pl . 



(9) 



Here Nf is a number of quark flavours, contributing to the as evolution. The 



onset of a new flavour q takes place at 



^4 



Ami, 



where, in general, m„ 7^ 
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Mg. Here are five different Ajg^(A^j) corresponding to Nj = 2,3,4,5,6 above tlie 
appropriate quark thresliolds. To sew up Og at tlie thresliolds an apparent matcliing 
condition for A's was used. Tlius at any ctsiQ"^) is determined by, say, as{M'^)^ 
wliich lias been used as a parameter for fits. 

We solved this equation at given numerically by rewriting it in the form x = 
f{x), where x = 47r/(/3oas)- It enabled us to obtain correct ag even at Q^/A^ ~ 1, 
where the expansion by powers of leading logarithms (see, e.g.,p^]) is no more valid. 
At Q^/A^ ^ 1 both methods yield the same result. (See Fig. 8) 



3.2 Preliminary fit results 

In order to check the consistency of the description of the data on 

R = a{e^e~ — > hadrons) / a{e^ e~ —>■ fJ.~^fJ.~) 

with the theoretical i2-ratio expression we have performed several fits using a sub-set 
of the total compiled data set of R which is applicable to the perturbative domain. 
This sub-set of the data is defined by the the following steps: 

1. The low energy region 2m,^ < ^/s < 2 GeV is completely excluded; 

2. All the data above 70 GeV are excluded, as we did not attempt to perform 
any fits of Z pole parameters Mz, Tz, s^r and Mto-p- This exclusion justified 
because these data slightly depend on as{Mz) and quark masses and thus 
their influence to the fits of these parameters would be rather negligable. 

3. Data located closer than 20 Breit-Wigner widths to the narrow hadronic 1 
resonances J/V'(1S), ij{2S), ^^(3770), T{IS), T(25), T(35), T(45), T(10860) 
and T (11020) are completely excluded; 



4. We also completely exclude SLAC-SPEAR-MARK-1 data ^ in the 2.6 GeV 
< < 7.8 GeV range. These data systematically lie ~15% above other 
experiments in the same interval. 

Some remarks on our treatment of the quark masses should be made here. We 
distinguished the 7* qq production threshold masses Mg and the QCD quark 
masses niq. The Mg's should be set actually equal to the masses of the lightest 
mesons, above which pair production threshold the onset of the new flavour q gives 
rise to the continuum cross section a{e~^e~ — > hadrons). The latters are just nig = 
fig/2, where fig is the energy scale at which the onset of the new flavour q in the 
evolution of Os takes place. We fitted only Mc and Mb, the other Mg and nig were 
fixed at their central PDG values. 
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The standard of the least squares method with weights as inverse squared 
total experimental errors (neglecting correlations in data) is used with the standard 
MINUIT[|lT5| package. 

After performing the exclusions 1) - 4) we fixed as{Mz), Mz, Tz, s'y^r, Mq 
{q = u, d, s, t), nig {q = u, d, s, c, b, t) at the values shown in the Table 1, leaving as 
free parameters just and Mf,. Another essential feature of our fit (I) was that 
we retained as a free parameters left and right boundaries of the interval [-^^1, ^/s2], 
3.0 < ^ < 3.670 GeV, 3.870 < ^ < 5.0 GeV, which was to be excluded from 
the data set in the process of /dof minimization itself.u Thus, the number of 
degrees of freedom was also variable during this fit. Such a /dof minimization 
procedure cancelled possible arbitrarities in the exclusions of broad resonances in 
the region above cc threshold. 

In the fit (II) we released as{Mz), fixing at the same time the excluded region 
^/s = 3.09 4.44 GeV. Mc, Mf, were fixed at their fitted values obtained in the fit 
(I). Fit (II) results in too high value as{Mz) = 0.128 ± 0.032. 

Preliminary fit parameter settings and results are enlisted in the Table 1. The 
reduced data set is shown on the Fig. 7. The references used for our fit and the list 
of excluded ^/s regions containing hadronic resonances are quoted in the Table 2. 



4 Summary. Assessed Data Compilations. 

In summary: 

• We have created two complementary computerized "raw" data compilations 
on a and R with data presented as in the original publications in all cases 
except the low energy subsample. In the low energy region, where there are 
no direct measurements of the total cross section, we obtain estimates of the 
total cross section either as the sum of exclusive channels or as the sum of the 
two-body exclusive channels with the data on e"*" e~ — >> Shadrons. 

• On the basis of above two data sets we have created compilations of data on 
^sd (corrected to the level of Fig 1) and B^"-"^^ (corrected to the level of Fig. 2) 
with one-to-one correspondence between the data points. Where necessary 

■^Narrow family resonances between 3 and 4 GeV were excluded by default. Somewhat controversial 
is the question, whether the gaps between t/i's can be really treated as the continuum regions. R ratio 
demonstrates no step up until = 3.9 — 4.0 GeV, just at the left knee of the broad ?/'(4040) resonance, 
where the continuum approximation of R does not work. Nevertheless, this step results in the fitted 
c-quark mass Mc — 1.971 (see fit (III) in the Table 1). 

In pQCD fits, performed in Refs. |l^, |^ the interval 3.0 < < 4.0 GeV with the resonances excluded 
was treated as a continuum. 
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the data have been rescaled to the standard style of implementing the pure 
QED radiative corrections to the initial state and to photonic propagator to 
produce a data set which is suitable for tests of the parton model with pQCD 
corrections for H^'"'^^ , and to be able to obtain more reliable estimates for 

• From the total data compilation on nf'"-^^ we have defined a "continuum" 
data compilation sub-set which can be used in conjunction with other data in 
the refinements of the as{Q'^) evolution and possibly in the global fits of the 
Standard Model. 

All data files are accessible by: 



http: / /wwwppds.ihep.su:8001 /comphp.html 



and will be accessible from the PDG site soon. 
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5 Appendix. Running aqED and Aag^^ 

As mentioned above, the authors of the original papers have published data in 
different forms which required the calculation of several types of factors, containing 
ctQEnis), to rescale the data to, say, "bare" i?-ratios. agEnis) can be expressed in 
the form (see: §, 0, §, M) 



aQEois) ^ ^^^^ 



l-Aa(s) ' 

Aa{s) = Aa^'^'^is) + Aa'^P(s) , (10) 
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where Aa'^"'^(s), Aa''^P(s) are hadronic and leptonic contributions to the QED vac- 
uum polarization, respectively. 

Aa''^P(s) is well determined in perturbative QED as a sum of leptonic loop 
contributions, 



a(0) 



In ^ - - + O ^ 



mf 



[ID 



The situation with Aa^'"^(s) is more complicated due to an essentially non- 
perturbative character of the strong interaction at low energy scales. Using the 
unitarity condition and the analyticity of the scattering amplitudes, one can express 
Aa^'^'^{s) in the form of a subtracted dispersion relation (see, e.g., ||5|, 0]) 



Aa^"'^(s) 



a{0)s R{s')ds' 
Svr J 4ml ^'i^' — s — iO) 



(12) 



where R{s) is the "bare" hadronic i?-ratio. 

This relation allows to effectively combine for its evaluation the pQCD /2-ratios 
in the continuum y/s intervals and experimental R data in the non-perturbative 
ones. 

We have evaluated the dispersion integral as 



a{0)s 



^ Rf'^ll{s')ds' 



Ami 



+ 



'Scut S'{s' — 



iO) 



(13) 



where the first integral was calculated numerically as a trapezoidal sum over the 
weighted average of experimental R points in the range 2m,r < < Jscut = 19.5 



GevH 

Our numerical evaluation procedure is in general similar to the one applied in 
Ref. 1^, exceptthat we performed the trapezoidal integration over all resonances, 
except (/.(1020).i 



^Choosing Scut we must satisfy at least the following two requirements: 

1) Scut should lie well above all the hadronic resonances, i.e. 1-loop perturbative QED can be used to 
calculate hadronic vaquum polarization at s > Scut] 

2) there should be enough densely distributed experimental points at s < Scut for the trapezoidal 
numerical evaluation of the dispersion integral over 2m^ < ^/s < y/s^. 

As there are few experimental points in the interval 13 < y/s < 30 GcV, s/Scut = 19.5 GeV appears to 
be a compromise between these requirements. 

^As the calculation of R'^^re from the original a and R data (see subsection 1.1) in turn requires the 
evaluation of Aa'""'(s), we applied the following simple iterative procedure. 



11 



We obtain here Aa'^'"^{M'^) = 0.0271 ± 0.0004(exp.), being consistent with the 
latest known to us results 0.027382 ± 0.000197 and 0.027612 ± 0.000220, obtained 
with two methodically different low energy data sets in Ref. pO[ |. 

The behaviour of Aa^^P, Aa['^^, q^^. and Aa'^^f . is depicted in Fig. 9. 



Taking as a zero approximation i?(s), obtained from the original data using 1-loop perturbative 
<3qed{s), we evaluate Aa'""*^^^ (s). Using the latter, we compute the next approximation of R{s) from the 
original data, the evaluate the integral again to obtain Aa'^'"''^) (s), and so on. This process converges 
well even in resonance regions after 3 iterations. We restricted ourselves to 5 iterations. (For such a 
number the finite machine precision is likely to result in the error, much less than the error of the method 
itself). 
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Table 1. Fit results in the defined "continuum" regions. Adjustable parameters are shown in 
bold. 



Parameter 


I 


II 


III 




U.iioi 






iviZ 


Q1 1S7 


Q1 1S7 


91 187 


J- z 


2. 1!) i i 


2.49-1 i 


2. 19 11 


a 


n 9'^! 1 7 




n 2*^1 1 7 


M, 




1 7A 


1 7A 


M 


n 140 


n i4n 


n 140 




n i/in 




c\^AC\ 


M 






n 4Q9 


M 




i.O 


x.y 1 ±u^± ) 


M, 


<^ O'i -1- 1 1 fi 




R (\ A- ^ A 
o.u zn 


rriu 


U.UUo 


U.UUo 


U.UUo 




n nnf; 

U.UUD 


n nnf; 

U.UUD 


n nnf; 

U.UUD 


iris 


n 1 on 


n 1 on 


n 1 on 
U.i/U 


rric 


1 


1 


1 o 
i.z 


Tin, I. 


4.2 


4.2 


4.2 




171.3 


171.3 


171.3 




'2m —20 


2?77 — 2 


Im -^2 




3.093-=-3.113 




3.093H-3.113 








3.684H-3.688 








3.670H-3.870 


Excluded 


3.175(247) ^ 


3.09 H-4.44 


4.000H-4.400 




4.319(105) 






Vs 


9.450^9.470 


9.450H-9.470 


9.450-=-9.470 




10.000h-10.025 


10.000-=-10.025 


10.000-=-10.025 


intervals 


10.34 H-10.37 


10.34 H-10.37 


10.34 H-10.37 




10.52 H-10.64 


10.52 H-10.64 


10.52 H-10.64 




10.75 H-10.97 


10.75 H-10.97 


10.75 H-10.97 




11.00 -=-11.20 


11.00 -=-11.20 


11.00 -=-11.20 




70 ^ 188.7 


70 ^ 188.7 


70 ^ 188.7 


xVdof 


0.690 


0.665 


0.822 
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Table 2. References of preliminar 



■y candidate data to form the continuum 



domain. 



Ref. No. 

Ill 
|69| 
1701 
|57| 
|33| 
|59| 
|27| 

1301 
|49| 
|39| 

m 

|46| 
|56| 
|37| 
|38| 
|44| 
1501 
|36| 
1401 
|35| 
|43| 
|52| 
|24| 
|31| 
|47| 
|26| 
|34| 
|55| 
|32| 
|51| 
|28| 
|25| 
1291 
|66| 
|67| 
|45| 



a type 


AT" 

points 




^max 


R 


85 


2.00 


4.80 


cr 


1 


2.23 


2.23 


cr 


11 


2.40 


5.00 


R 


6 


2.60 


5.00 


R 


2 


3.598 


3.886 


R 


33 


3.6025 


5.1950 


R 


27 


O O TO 

3.878 


4.496 


R 


15 


5.00 


7.40 


R 


31 


7.30 


10.29 


R 


3 


7.440 


9.415 


cr 


13 


9.30 


9.48 


R 


1 


9.36 


9.36 


R 


1 


n on 

9.39 


A OA 

9.39 


R 


1 


9.51 


9.51 


R 


1 


10.04 


1 A A /I 

10.04 


R 


1 


1 A AO 

10.43 


1 A AO 

10.43 


R 


1 


10.49 


in Ad 

10.49 


R 
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Table 3a. 



ZFITTER rescaling results for the Z pole. 
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Table 3a (continued). ZFITTER rcscaling results for the Z pole. 
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Table 3b. ZFITTER rescaling results for the LEP II-III data. 





GeV 






a^'^P [mb] 


^^^^•(^6orn/^cn*) [mb] 


130 


00000 


0.82691E- 


01 


0.82628E-01 


0.84200E-07 


0.84136E-07 


130 


11999 


0.82294E- 


01 


0.82264E-01 


0.79700E-07 


0.79670E-07 


130 


19999 


0.82031E- 


01 


0.82022E-01 


0.82100E-07 


0.82091E-07 


130 


19999 


0.82031E- 


01 


0.82022E-01 


0.71600E-07 


0.71592E-07 


136 


08000 


0.66194E- 


01 


0.67219E-01 


0.66700E-07 


0.67732E-07 


136 


10000 


0.66150E- 


01 


0.67177E-01 


0.66600E-07 


0.67634E-07 


136 


19999 


0.65930E- 


01 


0.66967E-01 


0.65100E-07 


0.66124E-07 


136 


19999 


0.65930E- 


01 


0.66967E-01 


0.58800E-07 


0.59725E-07 


161 


30000 


0.34218E- 


01 


0.35777E-01 


0.37300E-07 


0.38999E-07 


161 


30000 


0.34218E- 


01 


0.35777E-01 


0.40900E-07 


0.42763E-07 


161 


30000 


0.34218E- 


01 


0.35777E-01 


0.29940E-07 


0.31304E-07 


172 


10000 


0.28336E- 


01 


0.29787E-01 


0.30300E-07 


0.31852E-07 


172 


10000 


0.28336E- 


01 


0.29787E-01 


0.26400E-07 


0.27752E-07 


172 


30000 


0.28244E- 


01 


0.29693E-01 


0.28200E-07 


0.29646E-07 


182 


69000 


0.24017E- 


01 


0.25342E-01 


0.23700E-07 


0.25007E-07 


182 


69999 


0.24014E- 


01 


0.25338E-01 


0.21710E-07 


0.22907E-07 


182 


69999 


0.24014E- 


01 


0.25338E-01 


0.24700E-07 


0.26062E-07 


188 


63000 


0.22095E- 


01 


0.23352E-01 


0.22100E-07 


0.23358E-07 


188 


69999 


0.22073E- 


01 


0.23330E-01 


0.23100E-07 


0.24415E-07 



17 



References 



G. Bonneau and F. Martin, Nucl. Phys. B 27 (1971) 381. 

D. E. Groom et ai, [Particle Data Group Collaboration], Eur. Phys. J. 
C 15 (2000) 1. 

N. Cabibbo, R. Gatto, Phys. Rev. Lett. 4 (1960) 313; 

N. Cabibbo and R. Gatto, Phys. Rev. 124 (1961) 1577. 

R. B. Nevzorov, A. V. Novikov and M. I. Vysotsky, JETP Lett. 60 (1994) 
399 ||hep-ph/9405390|] . 



H. Schopper, D. R. Morrison, V. V. Ezhela, Y. G. Stroganov, 
O. P. Yushchenko, V. Flaminio and M. R. Whalley, Berlin, Germany: 
Springer (1992) 332 p. (Landolt-Boernstein. New Series, 1-14)- 

S. Eidelman and F. Jegerlehner, Z. Phys. C 67 (1995) 585 ||hep- 
ph/9502298| 



H. Burkhardt and B. Pietrzyk, Phys. Lett. B 356 (1995) 398. 
M. L. Swartz, Phys. Rev. D 53 (1996) 5268 ||hep-ph/950924g| . 



R. Alemany, M. Davier and A. Hocker, Eur. Phys. J. C 2 (1998) 123 
||hep-ph/ 9703220 1. 

M. Davier, LAL 98-87 (1998), Invited talk given at the "5*^ Int. Work- 
shop on T Physics," Santander (Spain), Sep. 14-17, 1998. 

M. Davier and A. Hocker, Phys. Lett. B 435 (1998) 427 ||hep-ph/9805"i70|] . 

S. Bethke, J. Phys. G 26 (2000) R27 [|hep-ex/0004021 



I. Hinchliffe and A. V. Manohar, Ann. Rev. Nucl. Part. Sci. 50 (2000) 
643 ||hep-ph/0004T8^] . 



Z.G. Zhao, Int. J. Mod. Phys. A 15 (2000) 3739. 
M. Vehman, Nucl. Phys. B 123 (1977) 89. 

J. Jersak, E. Laermann and P. M. Zerwas, Phys. Lett. B 98 (1981) 363. 
J. Schwinger, "Particles, Sources And Fields." Addison- Wesley, 1970. 3v. 
R. Marshall, Z. Phys. C 43 (1989) 595. 



A. D. Martin and D. Zeppenfeld, Phys. Lett. B 345 (1995) 558 |hep^ 
ph/9411377|| . 



18 



[20] A. D. Martin, J. Outhwaite and M. G. Ryskin, Phys. Lett. B 492 (2000) 
69 ||hep-ph/ 0008078 . 



[21] K. G. Chetyrkin, J. H. Kuhn and A. Kwiatkowski, Phys. Kept. 277 (1996) 
189. 

[22] K. G. Chetyrkin, R. V. Harlander and J. H. Kuhn, Nucl. Phys. B 586 
(2000) 56 ||hep-ph/ 00051391 . 



[23] D. Bardin, P. Christova, M. Jack, L. Kahnovskaya, A. Olchevski, S. Rie- 
mann and T. Riemann, Comput. Phys. Commun. 133 (2001) 229 hep- 
ph/9908433| 



bare 



R-ratio data with proper radiative corrections. 

[24] H. J. Behrend et al, [CELLO Collaboration], Phys. Lett. B 183 (1987) 
400. 

[25] H. Yoshida et al, [Venus Collaboration], Phys. Lett. B 198 (1987) 570. 

[26] E. Fernandez et al, [SLAC-PEP-006, MAC collaboration], Phys. Rev. D 
31 (1985) 1537. 

[27] A. Osterheld et al, [SLAC-SP-030, CRYSTAL BALL Collaboration], 
SLAC-PUB-4160. 

[28] 1. Adachi et al, [TOPAZ Collaboration], Phys. Lett. B 234 (1990) 525. 

[29] T. Kumita et al, [AMY Collaboration], Phys. Rev. D 42 (1990) 1339. 

[30] C. Edwards et al, [SLAC-SP-030, Crystall Ball Collaboration], SLAC- 
PUB-5160. 

[31] W. Braunschweig et al, [TASSO Collaboration], Z. Phys. C 47 (1990) 
187. 

[32] H. J. Behrend et al, [CELLO Collaboration], Phys. Lett. B 144 (1984) 
297. 

[33] P. A. Rapidis et al, [SLAC-SP-026, SMAG-LGW Collaboration], Phys. 
Rev. Lett. 39 (1977) 526. 

[34] D. P. Barber et al, [DESY-PETRA, MARK- J Collaboration], Phys. Lett. 
B 85 (1979) 463. 

[35] D. P. Barber et al, [DESY-PETRA, MARK- J Collaboration], Phys. 
Rept. 63 (1980) 337. 



19 



R. Brandelik et ai, [TASSO Collaboration], Z. Phys. C 4 (1980) 87. 

P. Bock et ai, [DESY-Hamburg-Heidelberg- Munich Collaboration], Z. 
Phys. C 6 (1980) 125. 

E. Rice et ai, [CESR-CUSB Collaboration], Phys. Rev. Lett. 48 (1982) 
906. 

B. Niczyporuk et al, [DESY-DORIS, LENA Collaboration], Z. Phys. C 
15 (1982) 299. 

R. Brandelik et al, [TASSO Collaboration], Phys. Lett. B 113 (1982) 
499. 

J. Siegrist et al, Phys. Rev. D 26 (1982) 969. 

L. Criegee and G. Knies, [PLUTO Collaboration], Phys. Rept. 83 (1982) 
151. 

B. Naroska, [DESY-PETRA, JADE Collaboration], Phys. Rept. 148 
(1987) 67. 

R. Giles et al, [CESR-CLEO Collaboration], Phys. Rev. D 29 (1984) 
1285. 

K. Abe et al, [VENUS Collaboration], Phys. Lett. B 246 (1990) 297. 

Z. Jakubowski et al, [DESY-DORIS, Crystal Ball Collaboration], Z. 
Phys. C 40 (1988) 49. 

C. Von Zanthier et al, [SLAC-PEP-005, MARK-II Collaboration], Phys. 
Rev. D 43 (1991) 34. 

H. Albrecht et al, [ARGUS Collaboration], Z. Phys. C 54 (1992) 13. 

A. E. Bhnov et al, [NOVOSIBIRSK- VEPP-4, MD-1 Experiment], Z. 
Phys. C 70 (1996) 31. 

M. Althoff et al, [TASSO Collaboration], Z. Phys. C 22 (1984) 307. 
M. Althoff et al, [TASSO Collaboration], Phys. Lett. B 138 (1984) 441. 

B. Adeva et al, [DESY-PETRA, Mark- J Collaboration], Phys. Rev. D 
34 (1986) 681. 

R. H. Schindler, [SLAC-SP-029, MARK II Collaboration], SLAC-R-0219. 

R. Ammar et al, [CLEO Collaboration], Phys. Rev. D 57 (1998) 1350 
||hep-ex/97070T^ . 



20 



D. P. Barber et al, [DESY-PETRA, MARK- J Collaboration], Phys. Lett. 
B 108 (1982) 63. 

H. Albrecht et al, [DESY-DORIS, DASP Collaboration], Phys. Lett. B 
116 (1982) 383. 

J. Z. Bai et al, [BES Collaboration], Phys. Rev. Lett. 84 (2000) 594 
||hep-ex/ 990801^ . 

J. Z. Bai et al, [BES Collaboration], |hep-ex/0102003 . 



R-ratio data, radiatively corrected according to Bon- 
neau/Martin. 

R. Brandelik et al, [DASP Collaboration], Phys. Lett. B 76 (1978) 361. 
cr data with proper radiative corrections, d^^. 

C. Berger et al, [PLUTO Collaboration], Z. Phys. C 1 (1979) 343. 

D. Andrews et al, [CESR-CLEO Collaboration], Phys. Rev. Lett. 50 
(1983) 807. 

R. Giles et al, [CESR-CLEO Collaboration], Phys. Rev. Lett. 50 (1983) 
877. 

J. K. Bienlein et al, Phys. Lett. B 78 (1978) 360. 

A. Boyarski et al, Phys. Rev. Lett. 34 (1975) 1357. 

J. Z. Bai et al, [BES Collaboration], Phys. Lett. B 355 (1995) 374. 

K. Abe et al, [TOPAZ Collaboration], Phys. Lett. B 304 (1993) 373. 

K. Miyabayashi et al, [TOPAZ Collaboration], Phys. Lett. B 347 (1995) 
171. 

cr data, radiatively corrected according to Bonneau/Martin. O"^. 

G. Cosme et al, [ORSAY-ACO Collaboration], Phys. Lett. B 48 (1974) 
155. 

B. Bartoh et a;.,[FRAS-ADONE, MEA Collaboration], Phys. Lett. B 58 
(1975) 478. 

J. E. Augustin et al, [SLAC-SP-017, SMAG Collaboration], Phys. Rev. 
Lett. 34 (1975) 764. 



a- data, from LEPI, LEPII-III. (jf. 



21 



[71] P. D. Acton et al, [OPAL Collaboration], Z. Phys. C 58 (1993) 219. 

[72] D. Buskulic et al, [ALEPH Collaboration], Z. Phys. C 60 (1993) 71. 

[73] P. Abreu et al, [DELPHI Collaboration], Nucl. Phys. B 417 (1994) 3. 

[74] P. Abreu et al, [DELPHI Collaboration], Nucl. Phys. B 418 (1994) 403. 

[75] M. Acciarri et al, [L3 Collaboration], Z. Phys. C 62 (1994) 551. 

[76] G. Alexander et al, [OPAL Collaboration], Z. Phys. C 52 (1991) 175. 

[77] G. Abbiendi et al, [OPAL Collaboration], Eur. Phys. J. C 6 (1999) 1 
||hep-ex/9808023|| . 

[78] P. Abreu et al, [DELPHI Collaboration], Eur. Phys. J. C 11 (1999) 383. 

[79] R. Barate et al, [ALEPH Collaboration], Eur. Phys. J. C 12 (2000) 183 
[|hep-ex/9904011|| . 

[80] G. Abbiendi et al, [OPAL Collaboration], Eur. Phys. J. C 13 (2000) 553 
[|hep-ex/9908008|] . 



^1] M. Acciarri et al, [L3 Collaboration], Phys. Lett. B 479 (2000) 101 |hep^ 
ex/000203l . 

cr data below 1 GeV. (Jg^. 

^2] I. A. Koop et al, [OLYA Collaboration], Budker INP-79-67 (1979) 

^3] A. Cordier et al, [ORSAY-ACO, DM1 Collaboration], Nucl. Phys. B 172 

(1980) 13. 

M] I. B. Vasserman et al, [VEPP-2, spark chamber system] Yad. Fiz. 33 

(1981) 709. 

^5] G. V. Anikin et al, [VEPP-2M, CMD Collaboration], IYF-83-85 Con- 
tributed to Int. Symp. on Lepton and Photon Interactions at High Ener- 
gies, Ithaca, N.Y., Aug 4-9, 1983. 

^6] L. M. Kurdadze et al, [VEPP-2M, OLYA Collaboration], Sov. J. Nucl. 
Phys. 40 (1984) 286. 

^7] V. L. Auslender et al., [VEPP-2, Spark chamber system], Yad. Fiz. 9 
(1969) 114. 

M. Bernardini et al, [FRASCATLADONE, BCF Collaboration], Phys. 
Lett. B 53 (1974) 384. 

[89] J. E. Augustin et al, [ORSAY-ACO], Phys. Lett. B 28 (1969) 513. 



22 



[90] D. Benaksas et al, [ORSA-ACO], Phys. Lett. B 42 (1972) 507. 

[91] P. M. Ivanov et al, [VEPP-2M, OLYA Collaboration], Phys. Lett. B 107 
(1981) 297. 

[92] D. Benaksas et al, [ORSA-ACO], Phys. Lett. B 39 (1972) 289. 

[93] M. N. Achasov et al, [VEPP-2M, SND Collaboration], |hep-ex/98090T3 . 

[94] M. N. Achasov et al, [VEPP-2M, SND Collaboration], Phys. Lett. B 462 
(1999) 365 ||hep-ex/991000l|] . 

[95] R. R. Akhmetshin et al, [VEPP-2M, CMD-2 Collaboration], |Ei^ 
ex/9904027| . 

References used for exclusive channels summation 
at energies 1.4 < Ecm < 2 GeV. cTq'^. 

[96] J. E. Augustin et al, [DM2 Collaboration], LAL-83-21 Contributed to 
Int. Europhysics Conf. on High Energy Physics, Brighton, England, Jul 
20-27, 1983. 

[97] D. Bisello et al, [DM2 Collaboration], Z. Phys. C 39 (1988) 13. 

[98] D. Bisello et al, [DM2 Collaboration], Phys. Lett. B 220 (1989) 321. 

[99] F. Mane, D. Bisello, J. C. Bizot, J. Buon, A. Cordier and B. Delcourt, 
[ORSA-DCI, DM1 Collaboration], Phys. Lett. B 99 (1981) 261. 

[100] I. A. Koop et al, [OLYA Collaboration], Budker INP-79-67 (1979) 

[101] L. M. Kurdadze et al, [OLYA Collaboration], Sov. J. Nucl. Phys. 40 
(1984) 286. 

[102] M. Bernardini et al, [ADONE, BCF Collaboration], Phys. Lett. B 46 
(1973) 261. 

[103] M. Grilli et al, [ADONE, MU-PI Collaboration], Nuovo Cim. ISA 
(1973) 593. 

[104] C. Bacci et al, [ADONE, 2GAMMA Collaboration], Phys. Lett. B 44 
(1973) 533. 

[105] C. Bacci et al, [ADONE, 2GAMMA Collaboration], Phys. Lett. B 58 
(1975) 481. 

[106] C. Bacci et al, [ADONE, 2GAMMA Collaboration], Phys. Lett. B 86 
(1979) 234. 



23 



[107] B. Esposito et al, [ADONE, MEA Collaboration], "Hadronic Cross- 
Section In E+ E- Annihilation From 1.45-Gev To 1. 80-Gev," Lett. Nuovo 
Cim. 30 (1981) 65. 

[108] B. Esposito et al, [ADONE, MEA Collaboration], "Multi - Hadron Pro- 
duction In E+ E- Annihilation At 1.45-Gev - 1. 61-Gev Cm Energy," Lett. 
Nuovo Cim. 25 (1979) 5. 

[109] B. Bartoli et al, [ADONE, BOSON Collaboration], Phys. Rev. D 6 
(1972) 2374 

[110] C. Bacci et al, [ADONE, 2GAMMA Collaboration], Phys. Lett. B 38 
(1972) 551 

[111] F. Ceradini et al, [ADONE, MU-PI Collaboration], Phys. Lett. B 42 
(1972) 501 

[112] S. I. Dolinsky et al, [VEPP-2M, ND Collaboration], Phys. Kept. 202 
(1991) 99. 

[113] jittp: / / wwwppds.ihep.su:8001/ ppds.html 
[114] [http: / / durpdg.dur.ac.uk/hepdata/l 

[115] F. James and M. Roos, Comput. Phys. Commun. 10 (1975) 343. 



Received 16 May, 2001 



24 



Figure Captions 

Fig. 3: World data on the total cross section of e~^e^ — > hadrons as pre- 
sented in the original publications. 

Fig. 4: World data on the experiment-to-theory ratio R = a{e^e~ — >■ 
hadrons) experimental/ (^{e+e~ fi'^fi~)QED simple pole as presented in the origi- 
nal publications. 

Fig. 5: World data on the total cross section of e"'"e~ hadrons. Curves 
are an educative guide. Solid curve is the cross section prediction obtained in 
the three-loop QCD approximation with the effect of non-zero quark masses 
taken into account. 



Fig. 6: i?-ratio. Data set is the same as on the Fig. 5. Solid curve is the 
i?-ratio prediction obtained in the three-loop QCD approximation with the 
effect of non-zero quark masses taken into account. Dashed curve is a "naive" 
quark parton model prediction for the i?-ratio 



Fig. 7: Data set used for the fits. 

Fig. 8: A comparison of i?-ratio obtained with as{s), evaluated by our 
numerical method (solid curve) and by the method described in [jl3| (dashed 
curve). 

Fig. 9: The behaviour of Aoqed- ^Q^iaJ^^ the quark contribution to 
Aoi(^QED) calculated in 1-loop QED. Aa^^^^^^ is a hadronic contribution to 
Aa(^QED) calculated using the dispersion integral over crhad{s), 2mT, < < 20 
GeV, and 1-loop QED above 20 GeV. 
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Figure captions 



Fig. 3: World data on the total cross section of e'^e hadrons as pre- 
sented in the original publications. 

Fig. 4: World data on the experiment-to-theory ratio R = a[e~^e~ 
hadrons)experimentai / o-{e~^ e~ fi'^ fJ'~)QED simple pole as presented in the orig- 
inal publications. 

Fig. 5: World data on the total cross section of e"'"e~ hadrons. Curves 
are an educative guide. Solid curve is the cross section prediction obtained in 
the three-loop QCD approximation with the effect of non-zero quark masses 
taken into account. 



Fig. 6: R-vatio. Data set is the same as on the Fig. 5. Solid curve is 
the R-vatio prediction obtained in the three-loop QCD approximation with 
the effect of non-zero quark masses taken into account. Dashed curve is a 
"naive" quark parton model prediction for the R-vatio 

Fig. 7: Data set used for the fits. 

Fig. 8: A comparison of R-vatio obtained with as{s), evaluated by our 
numerical method (solid curve) and by the method described in [13] (dashed 
curve). 

Fig. 9: The behaviour of AaqED- •^ctiad'^^ quark contribution to 

Aci{QED) calculated in 1-loop QED. Aaj^^^^^ is a hadronic contribution to 
Aa(^QED) calculated using the dispersion integral over (Thad{s), 2m^ < < 
20 GeV, and 1-loop QED above 20 GeV. 
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